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' Abstract. There are several methods to treat ensembles of random matrices 

in symmetric spaces, circular matrices, chiral matrices and others. Orthogonal 
F ^ | polynomials and the supersymmetry method are particular powerful techniques, 

rj^ ■ Here, we present a new approach to calculate averages over ratios of characteristic 

. polynomials. At first sight paradoxically, one can coin our approach "supersymmetry 

without supersymmetry" because we use structures from supersymmetry without 
' ' actually mapping onto superspaces. We address two kinds of integrals which cover 

O j. a wide range of applications for random matrix ensembles. For probability densities 
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factorizing in the eigenvalues we find determinantal structures in a unifying way. As 
a new application we derive an expression for the fc-point correlation function of an 
arbitrary rotation invariant probability density over the Hermitian matrices in the 
presence of an external field. 
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1. Introduction 

OV 

Random matrix theory has a wide range of applications [U El El E] resulting in a large 
number of different matrix ensembles. To name but a few examples, generic features of 
Hamilton operators are modeled by ensembles over the symmetric spaces [3 [21 [6] . In a 
relativistic setting, chiral (Laguerre) random matrix ensemble have to be used, involving 
one [31 [7] or two [H El [TUl HH [12] matrices. Another example is the statistics of the 
zero points of the Riemann zeta function [TSJ, Q31 [T51 dS] , derived by Dyson's circular 
ensembles [17] . 

Averages over ratios of characteristic polynomials play an important role in the 
investigation of random matrix ensembles. The matrix Green function can be generated 
by one characteristic polynomial in the denominator and one in the numerator [T8l IB]. 
For the calculation of the free energy, one may use the replica trick [TH]. The 
moments of the Riemann ^-function are also of interest for number theorists [14, 15J. 
Mathematicians are interested in averages over ratios of characteristic polynomials 
because of the connection to Weyl's character formula [201 EI]- In models for Quantum 
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Chromodynamics (QCD) [3], [9] and in the analysis of the sign problem [22], one employs 
mean values of characteristic polynomials. 

To calculate such mean values, commonly two techniques are used, the method of 
orthogonal polynomials [221 EU H] and the supersymmetry method [251 El El ESI ETJ 
1281 l29j [30] . In both methods determinantal structures appear for rotation invariant 
ensembles over the Hermitian matrices [311 E21 [33J EH [35]. The probability densities of 
such ensembles have to factorize in the eigenvalues of the matrices. The determinantal 
structures extremely simplify the calculation since all fc-point functions are determined 
by one-point and two-point correlations. Also for other ensembles such as the non- 
Hermitian ensembles [36j l2l] . the chiral ensembles [HOE] and the circular ensembles 
[3T1 EB] , determinantal structures were found. 

In the orthogonal polynomial method as well as in the supersymmetry method 
every single ensembles has to be calculated in a particular way. Either one has to 
find the measure to construct the orthogonal polynomials or one has to identify the 
superspace corresponding to the ordinary integration domain. We consider two types of 
integrals related to mean values of ratios of characteristic polynomials. Determinantal 
structures stemming from supersymmetry, such as those found by Basor and Forrester 
[37] , yield determinantal structures of these integrals. Here, we establish the link to 
supersymmetry. To the best of our knowledge this connection has not been observed 
before. Our method is based on an algebraic manipulation of the characteristic 
polynomials and the Jacobian or the Berezinian resulting from changing integration 
variables. We neither use the Mehta-Mahoux theorem nor a mapping onto superspace. 
Both types of integrals cover a wide range of applications for unitarily rotation 
invariant random matrix ensembles and for ensembles over all three rotation groups 
or their Lie algebras. As a particular example, we consider an intermediate ensemble 
from arbitrary unitarily invariant ensembles over Hermitian matrices to a rotation 
invariant ensemble over one of the symmetric spaces. This generalizes known results 
Sni SU SSI S3l SSI SSI SZ] - For this example we use the supersymmetry method. 
Thereby, we demonstrate that our method works for calculations within superspace, 
too. 

In Sec. [21 we give an outline of our approach. We present a determinantal 
structure of Berezinians resulting from a diagonalization of symmetric supermatrices 
in Sec. EJ This is then applied to two types of integrals discussed in Sec. H] which yield 
determinantal structures. In Sec. El we present a series of ensembles whose mean values 
of characteristic polynomial ratios are special cases of one type of integral presented in 
Sec. [H In particular, we consider Hermitian matrices which we investigate in Sec. [6] 
as well. In Sec. [6] we, also, study examples for the other type of integral. The main 
focus is on the Hermitian matrix ensemble with arbitrary unitarily rotation invariant 
probability density in the presence of an external field. In the appendices, we perform 
some explicit calculations. 
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2. Sketch of the idea: "Supersymmetry without supersymmetry" 



As a guideline for the reader, we present here the main ideas of our approach for one 
particular ensemble. We choose k = diag (/tn, . . . , k 12 , . . . , « fc2 ) = diag k 2 ) 
in such a way that the integrals below are well defined. For many applications such 
as for Hermitian matrix ensembles one considers averages over ratios of characteristic 
polynomials 

k 

det(H — k^Ijv) 



Z k) 



P(H)\\ 



-d[H] . 



(2.1) 



11 det(H 
3 =i 

Here, Ijv is the N x N unit matrix. The probability density P is rotation invariant 
and factorizes in the eigenvalues of the matrix H. We diagonalize H in its eigenvalues 
Ei, ... , En. The Jacobian is the second power of the Vandermonde determinant A^(E). 
We expand one of the Vandermonde determinants and have up to a constant c 



Z(k) 



N 

n 

a=l 



p{E a )E^ n 



6=1 



A Af (E)d[E] . 



(2.2) 



We pursue an idea similar to the one by Basor and Forrester j37j. We supplement the 
factor 

N k 



— Kbl 



0=1 6=1 



by 



BergU^) 



A fc (Ki)A fc (K2) 



(2.3) 



(2.4) 



n (« a i - 

a,6=l 

Both factors together are up to a sign 



Ber 



(2) 

k/k+N 



[ki, k 2 ,E) 



A k (Ki)A k+N (K 2 ,E) 



' k k N 

n (^ai - K b2 ) n n («oi 

a,6=l a=l 6=1 



(2.5) 



^6 



The authors in Ref. [37] have shown that this function has for all iV G N a determinantal 
structure mixing terms of the Vandermonde determinant and the Cauchy determinant, 



Ber 



(2) / 
k/k+N^ 





r i 


1 


k 2 , E) = ± det 


«al — K 62 


Koi — Eb 




K 62 


El 1 



(2.6) 



The insight crucial for this work and not contained in Ref. [37J is the intimate connection 



of Eq. (J2SD to superspace: Ber^^ is the Jacobian or Berezinian for the diagonalization 



of symmetric matrices or supermatrices, respectively, as shown in Refs. jH] . We will 
prove Eq. (12.61) in a new way and obtain also an interesting intermediate result not given 
in Ref. [37], see Sec. |3j 
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We now proceed with the evaluation of Z(k). We shift the eigenvalue integrals into 
the determinant and obtain 



Z(k) 





r i 


F b {K al ) 


det 


Kal — K b2 




K b2 


M ab 



(2.7) 



The symmetric matrix M ab comprises the moments of the probability density P and the 
functions F b are the Cauchy transform of those moments. At this point we have a choice 
for how to proceed further. For instance we can reorder the monomials in the entries 
of the determinant to orthogonal polynomials with respect to P. Then, M ab becomes 
diagonal and F b are the Cauchy transforms of the orthogonal polynomials. Thus we 
arrive at the well known result, see Ref. (32]. On the other hand, we can choose an 
arbitrary set of linearly independent polynomials. Then, we use the important property 
of the determinant 



det 



.4 
C 



B 
D 



det D det[A - BD^C] 



(2.{ 



for arbitrary matrices A, B and C and an invertible quadratic matrix D. This finally 
yields 

n 



Z(k) 



det 



«al — K b2 



- ]T F m {n al )M-^n n b - 1 



m,n=l 



det K(k u i, K b2 ) 



(2.9) 



We obtain the correct result without the Mehta-Mahoux theorem for an arbitrary 
choice of polynomials. The orthogonal polynomials are not the tool to identify the 
determinantal structures. They are a result of the calculation. 

In the next sections, we extend this sketch to a careful discussion for a large class 
of integrals. We will see that determinantal structures derived in many different fields 
of random matrix theory have a common origin. 



3. Determinantal structure of Berezinians 



In Sec. I3.1[ we investigate the determinantal structure of the Berezinians resulting 
from Hermitian supermatrices. These Berezinians are crucial for the calculations in the 
ensuing sections. For the sake of completeness, we give the determinantal structure 
according to the supergroup UOSp (p/q) in Sec. 13.21 

3.1. Berezinians related to the supergroup U(p/g) 

As we have seen in Sec. [2, Berezinians resulting from diagonalization of supermatrices 
play a role of paramount importance for our method. Although we do not use 
any integral in superspace we find those Berezinians in the ratios of characteristic 
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polynomials times the Vandermonde determinant. The crucial step is here to understand 
that those Berezinians have always a determinantal structure. 

For the Vandermonde determinant the determinantal structure 

A fe (K)= J] (Ka-« 6 ) = (-l)*t*- 1 ^det[«r 1 ] x < a «*. (3- 1 ) 

l<a<6<fc 

has been known for a long time [50j. Moreover, the square root of the Berezinian 
resulting from a diagonalization of the supersymmetric analog of a (k + k) x (k + k) 
Hermitian matrix is known |48j, up to a sign, to be equal to the Cauchy determinant 



n («al - Kbl) («o2 - ^fe2) 
l<a<fe<fc 

n (/« a i - ^62) 

Ka,£><fc 



. 1)fc (fc-l)/2 det 



1 



«al K b2\i <ab<k 



(3.2) 



The last equality sign is due to Cauchy's lemma [51J. The upper index "(2)" refers to 
the Dyson index e {1, 2, 4} and indicates that this Berezinian is related to Hermitian 
supermatrices. 

The next step is to generalize these structures to an arbitrary number of bosonic 
eigenvalues K\ = diag («n, . . . , k p i) and of fermionic eigenvalues K2 = diag (k-u, • • • , Kqz)- 
In Appendix A.l we derive the determinantal structures of the Berezinians, see Ref. 

EI - K&l) EI ( K a2 ~ K b2 ) 

l<a<b<p l<a<b<q 

p q 

n n - ^2) 

0=1 6=1 

for arbitrary p and q. Equation ( 13. 2 p is the case p = q = k. Under the condition p < q, 
we obtain 

p-q q-p 
^al ^62 



Ber^(^) 



(3.3) 



Ber( 2 / ) o ( K ) = (-l)^- 1 )/ 2+ (' ?+1 ) p det 



Kb2 J !< a <P 
l<b<q 



l<a<q— p 
l<b<q 



(3.4) 



Since the left hand side is up to (— l) pq symmetric under exchanging the bosonic 
eigenvalues with the fermionic ones, the condition p < q is not a restriction. This 
result is similar to Eq. (12. 6p . In Sec. 16.21 and Appendix E we show that it is useful for 
some calculations. 

The left hand side of Eq. (13. 4p is translation invariant k — > K+el p+q with a constant 
e. Thus, we may shift the expressions on the right hand side by e, 



Ber^(/c) = (-lyh-VM**- 1 )? det 



{i^ + ey- 1 } 



l<a<p 
l<b<q 



l<a<q— p 
l<b<q 



(3.5) 
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We expand the entries in the lower (g — p) x q block in e. We notice that all rows 
together are a linearly independent set of polynomials from order zero to order q—p — 1. 
As the determinant is skew symmetric, it yields 



Ber^(«) = (-1)9(<?- 1 )/ 2 +(«+ 1 )p det 



K>al — K b2 
{<2 1 } 



l<a<p 
l<b<q 



l<a<q— p 
l<b<q 



(3.6) 



Since e is arbitrary we take the limit for e to infinity and obtain the final result 

1 



Berfl(K) = (-lyte-W+it+VP det 



K>al — K b2 J l<a<P 
l<b<q 



<a<q—p 
l<b<q 



(3.7) 



which is identical to the result of Basor and Forrester [37] . Indeed, Eq. (13. 7p does not 
exhibit the nice symmetry between the bosonic and fermionic eigenvalues as in Eqs. ( 13. 2 j) 
and IEE51). 



3.2. Berezinians related to the supergroup UOSp(p/g) 

As for the Vandermonde determinant itself, the determinantal structure for the forth 
power thereof is also well known [50] . 



l<a<b<k 

det 



(a - Ik 



a-2 



l<a<2fc 
Kb<fc 



(3.8) 



Recently [29], it was shown that the Berezinian corresponding to the supergroup 
UOSp (2k /k) has a determinantal structure, too, 



Ber 



(i) 

2k/k 



k) = Ber 



(4) 
k/2k 



EI - «6l) 11 ( K a2 ~ ^2)" 

l<a<6<2fc l<a<fe<fc 

2k k 

II II ( K al ~ Kb2f 
0=1 6=1 



det 



«ol — K 62 



,« ol - K b2) 



l<a<2k 
Kb<k 



(3.9) 



Here, we define k = diag (k 2 , 



In Appendix A. 2 we derive the analog of Eq. (13. 3p for the Berezinian 

Ber ( V (k) = Ber (4 ; = 



p q 

n n (^ai - K b2 ) 

a=l 6=1 



(3.10) 
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see Ref. [53]. Here, we have to distinguish between p < 2q and p > 2q. For the first 
case, we obtain 



Ber^,(«) 



Ber$(*) 



-l?det 



K P-2q 2q-p 
„a-l 



1,2 



r) K P ~ 2q K 2q ~ p 

U K al K b2 
dKb2 K al — K b2 

(«-l)< 2 " 2 } 1 



l<a<p 
l<b<q 



<a<2q— p 
l<b<q 



(3.11) 



and, for the second one, we get 



Ber«(K) 



-l)p(p-i)/2+a de t 



(3.12) 



Kgl ^62 



Kgl ^62 
(«al - «62) 2 



l<a<p 
l<ft<c? 



Rt 1 } 



l<a<p 
l<6<p-2g 



We apply the same procedure as in Sec. 13. II and shift all elements by e. Taking the limit 
e — > oo, we find 



Ber^) 



BerSJ,(«) 



-l) p det 



1 



Kal ~ K b2 
{ K b2 1 



1 



(a-l)<~ 2 } 



l<a<p 
l<b<q 

l<a<2q-p 
l<b<q 



(3.13) 



Ber (4 ? (k) 



for p < 2q and 

Ber (1 ? (k) = Be, , 

p/? v 1 q/p 

(-l) p(p - 1)/2+ « det 



(3.14) 



Kal — Kb2 



[Kal - K h 2) 



l<a<p 
l<b<q 



<P 
l<6<p-2ij 



for p > 2q. We notice that the determinantal structure of the ordinary Jacobians 
which are powers of Vandermonde determinants mixes with the structure of the Cauchy 
determinant, as for the U (p/q) case. 



4. Main result 



We discuss two types of integrals which cover averages over ratios of characteristic 
polynomials for a large class of matrix ensembles. For both types we find determinantal 
structures. These types are integrals with a square root of a Berezinian f 1 3 . 3 p and with a 
Vandermonde determinant to the second power. They are studied in Sees. 14.11 and 14. 2\ 
respectively. 
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4-1. Integrals of square root-Berezinian type 



We consider the integral 

Ni N 2 



k 1 /k 2 V 



X 



3=1 



3=1 



£N 1 +N 2 

Ni fci N 2 k 2 

n n ( z ai - «6i) n n (z a2 - K b2 ) 

0=1 6=1 a=l 6=1 

Ni k 2 ki N 2 

n n ( z ai - K b2 ) n n (^1 - ^ 

a=l 6=1 a=l 6=1 



Ber !v!/7v 2 ( z ) rf w > 



(4.1) 



where the Zj are complex variables. The functions gj and fj and the variables k are 
chosen in such a way that the integral is convergent. The measure d[z] is the product of 
the differentials of the real and imaginary parts. The applications which we will give are 
particular cases of this integral, although these applications correspond to essentially 
different ensembles. Thus, we show a fundamental relation which yields determinantal 
structures. 

The crucial step is to extend the integral (14.11) by ^/r3er ^/ fe2 (ft) and to recognize 
that we obtain a new Berezinian 

iVi N 2 

' N *)f^ - I 
ki/k 2 



3=1 



(4.2) 



£N 1 +N 2 

where the new bosonic eigenvalues are z\ = diag (z±, Ki) and the new fermionic 
eigenvalues are z 2 = diag (z 2 , k 2 ). Now we use the determinantal structure of the square 
root Berezinian shown in Sec. 13.11 

Without loss of generality, we assume N 2 > N\. In Appendix B.2, we show the 
details of this calculation and only give the results here. The simplest case is k\ = k 2 = k. 
Then, the condition (Ni + ki) < (N 2 + k 2 ) is automatically fulfilled. The integral 
is then a quotient of two determinants times a constant 



where we define 



(_l)(" a +*)(i\r 3 +*-i)/2 det Mni/N2 



^k%(K) 



det [K^\n al ,K b2 )] 1 



<a,b<k 



(4.3) 



G^\n b2 ) 



<N 2 ), 



K a l) 



Z — K b2 



l<a<N 2 -N x 



-d[z] 



Ka<N 



d[z] 



l<b<N 2 



(4.4) 



(4.5) 
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Ni/Na 




l<a<N 2 -N 1 
l<b<N 2 



l<a<N 1 
l<b<N 2 



Kb2) 



F^\ Kal )M N ] /N G 



(N\ I N 2 ) / 



(4.6) 



(4.7) 



K al — K a2 

Since the entries K^ Nl ^ N2 ^ (K a i, ^2) are independent of the dimension k, we identify 

^W,^ ) = (=!^g^ (4,, 

det MjVi/ATa «ol - K b 2 

which is the case k = 1. The normalization constant follows from k = and is given by 
C Ni/ n 2 = Z$ /N2) = det M Nl/N2 . (4.9) 

This leads to the very compact result 



7 {Ni/N 2 ) 
J kjk 1 



(_l)fc(fc-l)/2 



/~ffc— 1 

°JVi/JV 2 



Ber i 2 ?, (re) 



det 



(N 1 /N 2 ) l 



7VJ 
Z l/1 



(4.10) 



Ka,Kt 



We recall that the functions gj and fj are arbitrary. This means that the fundamental 
structure of the ratios of the characteristic polynomials times the Berezinian fully 
generates the whole determinantal expression. 

The cases k\ < k 2 and k\ > k 2 cover all cases mentioned above. As in Ref. [8], 
we trace these cases back by introducing \k\ — k 2 \ dummy variables. These variables 
enlarge the {k\ + k 2 ) x (k\ + k 2 ) eigenvalue matrix of k to a {k + k) x {k + k) eigenvalue 
matrix, where k = max{fc l5 fc 2 }- Then, we use our result obtained above and remove 



these additional eigenvalues. The explicit calculations are given in Appendix B.2 and 



Appendix B.3 We obtain 



'k x ik 2 ' 



(_l)fei(fci-l)/2+(fe 2 -fei)iVi 



rihi—\ 
°JVi/JV 2 ' 



'Ber 



(2) 



X 



x det 




ki/k 2 

(iVi/iVa) 
1/1 



-1 ^iVa-ATi+l 7 {N 1 /N 2 ) 
K ^1/1 1 



l<a<fci 
l<6<fc 2 



K 0, K 62; 



^0 — K b2 



Kg— ¥00 ■ 



l<a<k 2 —k\ 
l<b<k 2 



(4.11) 



(_l)(fc2+2fc 1 )(fc 2 -l)/2+(fc 1 -fc 2 )(A r 2-iV 1 ) 
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x det 



^1/1 l K 61;^a2j 



.2 ± 



Kbl ~ K a2 

a-1 L ,N 1 ~N 2 + l 7 (N 1 /N 2 ) 



l<a<k 2 
l<b<k! 



'0 



'1/1 



l<a<fci— k 2 
no-tocs l<b<k! 



(4.12) 



for ki > k 2 . For \k% — k 2 \ = 1 the average in the last row is only over one characteristic 
polynomial, i.e. it is equal to Z^J^ ' \k\o) in Eq. (14.111) and Z±?^ N2 \Kbi) in Eq. (I4.12p . 
The limits in Eqs. (14. lip and (I4.12p are well defined, as a comparison with Eq. (14. ip 
shows, and can be calculated by writing the derivative as a contour integral around 
1/kq = 0. The limits are explicitly given as 



g 
dK 



a-1 .^-JVl+l/W/^), 



'1/1 



^0, K b2) 



K — K b2 



KQ—^OO 



"I) 



o-l+AT 2 , 



a - 1)!C. 



N1/N2 



K 



a-l+N 2 -N- 
b2 



^-i a M N l l/N G^\K b2 ) 



(4.13) 



and 







-1 .JVi-iVa+l y(N 1 /N 2 ) 
K ^1/1 



v«61, «0j 



/ K bl ~ K 



(-l) a+N *(a-l)\C Nl/N2 x 



X 



K 



a-l+Ni-N 2 
bl 



e( a + iVx - JV a - 1) - F^(«6i)M^ /JV3 ga 



where we define the matrices 



a-l+JV 2 -iVi 



d[z) 



l<b<7V 2 

{— 5jv 2 -JVi+l-a,6/ 

' l<fe<V 2 -Vi 

:)z a - 1+Nl - N2 d[z]e(a + Nt - N 2 - 1) 



l<fe<Vi 



(4.14) 



(4.15) 



. (4.16) 



The function is the Heaviside distribution for discrete numbers which means it is the 
integrated Kronecker-5 and, hence, unity at zero. 



4-2. Integrals of squared-Vandermonde type 
Now, we investigate integrals of the type 



h/h 



C N 



k 2 N l 2 N 

n n n n n(Ao2-^) 

n^-)^¥ ^ \^{z)?d[ Z] 

U U(^l-Zb)U n(Aal-4) 



(4.17) 



a=l 6=1 a=l fe=l 

The function g is, as in Sec. I4.1[ an arbitrary function with the only restriction that 
the integral above is convergent. Instead of one eigenvalue set as in the subsection 
above, we have now two eigenvalue sets n = diag (ku, . . . , k^i, K12, ■ ■ ■ , «fc 2 2) and 
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A = diag (An, . . . , A^i, A12, • • • , A^ 2 2). Because of these two sets, we have to extend 
the fraction by two square roots of Berezinians and find 



HN) 

' fel/*2 

h/h 



K, A) 



N 



Ufa) 



Beif i/k2+N (~z) 



We introduce 5 = diag («i, k 2 , z) and 5 = diag (Ai, A2, z*). 

To integrate Eq. ( I4.18p . we first discuss the case d = k 2 + N — k\ — l 2 + N — l\ > 0. 
Under the integral we have two determinants with N rows depending on one z a or one 
z*. The other rows are independent of any z a and z*. Thus, we use an integration 
theorem similar to Andreief 's [SI] which we derive in Appendix C.l In Appendix D we 
carry out the calculation and find 

/ -l)(b+fe 2 )(Ji+fci-l)/2 iV | det M (i 



Ber 



(2) 



(A) 



d[z\ 



(4.18) 



h/h 



K^f(n a i, K b2 ) 



x det 



where 



l<a<Z 2 
l<ft<fc 2 

l<a<fci 
l<6<fc 2 



M/0 ( K al; Afti) 



1/0 
1/0 



BergjA) 

^12^(Abl) A a 2 
K 22 \n a l, Aw" 



X 



(Kal 



Fd(«ol! 



^)(A 6 1 



l<a<« 2 
l<6<Zl 

l<a<fci 
l<fe<«i . 



-d 2 z , 



F d W (A 



61, 



M d 

Ad(A a2 ) 

K d (AC 62 ) 

A^f (k 6 2, A a2 ) 

K^(Xbl, A a 2) 

K 2 l( K al, K b2 ) 
K 2 f(n a l, Afcl) 



5f(^)2; a_1 
Aw - 2* 



d 2 z 



J K6<d 



Ka<d 



g{z)z a - 1 z* b - 1 d 2 z 



Ka,b<d 



Ns 1 ] 

[ K 62 ] l< a <d > 

- A d (A aa )M^ 1 K < ,(K6 2 ) , 

t — — r A d (A a2 )M d 1 F]* ) (A w ) , 

^W — ^a2 

F d (K al )M^K d (K b2 ) 

K-al — K> b 2 

Z ( ^{K al ,X hl ) - ¥ d {K al )m- d l ¥ { ;\x bl ) 

1/0 



(4.19) 



(4.20) 
(4.21) 
(4.22) 

(4.23) 

(4.24) 
(4.25) 
(4.26) 
(4.27) 

(4.28) 
(4.29) 
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With help of the particular cases (ki — k 2 — 1, l\ — l 2 — 0), {ki = k 2 = 0, l\ = l 2 = 1) 
(ki = li = 1, k 2 = l 2 = 0) and (fci = l\ — 0, k 2 — l 2 — 1), we identify 



^1/1 ( K al: K b2) 



III 

0/0 



a2j 



-"■22 [ K al,A b l) 



JVIdetMjv K al -K b2 

0/0 

1/1 



1 



^0/0 (A&1> ^02) 



N\ det M N 
1 



K 



(N+l) 
11 



(A a 2, K b2 ) 



iV!detMjv_i 1/0 
1 



^1/0 { K al, 



bh 



7 (JV) / \ \ 
^0/1 l^a2, «b2j • 



iV!detMjv + i 0/1 

The normalization constant is fixed by the case {k\ = k 2 = l\ = l 2 = 0) 



C 



N 



1 0/0 
0/0 



N\ det M 



N 



(4.30) 

(4.31) 
(4.32) 
(4.33) 

(4.34) 



Thus, we find 



h/h 



(-l\{h+k 2 )(h+ki-l)/2 N ] 



x det 



det^ 1 " 1 

Zq/i \^a2i K b2) 
0/1 



^0/0 (Afei? A a2 ) 
1/1 



l<a</ 2 
l<b<fe 2 



(d-1)! 

M/i (^alj ^62) 
0/0 

d!(/C a l - K62) [ l<a<fci 
!<£><fc 2 



d!(A 6 i - A a2 ) [ l<a</2 

l<b<h 

z[/t ( K <zi> A{,i) 



1/0 



(tZ+1)! 



l<a<fci 
l<b<h 



(4.35) 



Once more, we notice that the distribution g(z) is quite arbitrary and, thus, a large class 
of ensembles is covered. The result (I4.35P is equivalent to the one found by Bergere [21] 
with the method of biorthogonal polynomials. 

We derive the integral (Q7j) for the case that d = k 2 + N-k 1 = l 2 + N-l 1 >0 
is violated by the same method as used in Sec. 14.11 By extending the quotient of 
characteristic polynomials to the case discussed above, we apply the known result and 
take the limits with help of l'Hospital's rule. This procedure gives us expressions similar 
to Eqs. fl4TTT]) and (Q21) . 

For the particular case that d K = k\ — k 2 — N and d\ — 1% — l 2 — N are larger than 
zero, we find a much simpler expression 

(_l)(h+k 1 )(h+k 1 -l)/2+N(k 2 +l2+l) Af\ 

A) = 1 } x 



Ber^ 2 («)VBer^ 2 (A) 
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x det 



Vi2 



Ka 



l<a<Z 2 
Kb<L 



1<6<«1 



l<a<fci 
l<6<fc 2 



<fcl 
Kb<d K 



1/0 



l^ali A 



61, 



l<a<fci 
l<fe<«i 



(4.36) 



which is derived in Appendix D.2 We see that the whole integral is determined by 

^i/o ( K ai, Afci) and some algebraic combinations of the "Cauchy terms" l/(\bi — A a 2) 

170 

and l/(ft a i — ^62) arid of the "Vandermonde terms" A^f 1 and k^ 1 if the number 
of characteristic polynomials in the denominator exceeds a critical value. Here, we 
emphasize that Eq. (I4.36P has a simpler structure than Eqs. (14. lip . (I4.12p and (I4.35P 
for d K 7^ d\ since there are no limits to perform. Furthermore, the parameters d\ and 
d K are independent as long they are larger or equal than zero. 

Eqs. (I4.35P and (I4.36P are similar to the results found by Uvarov [551 ES]- He 
studied the transformation behavior of the orthogonal polynomials when the probability 
distributions differ in a rational function. This shows the connection between his 
approach and ours. 



5. Applications for integrals of squared— Vandermonde type 

In Sec. 15.11 we apply our method for integrals of squared- Vandermonde type to the 
example of Hermitian matrices. Since the method has a broad field of applications, we 
give a list of matrix ensembles in Sec. 15.21 



5.1. Hermitian matrix ensemble 

We first consider rotation invariant ensembles over the N x N Hermitian matrices 
Herm(V). The averages 

J! det(if - K j2 t N ) 
4f/U-)= J P (H)^ d[H] (5.1) 

Herm(V) FJ det(iJ - Kjil N ) 

3=1 

are of considerable interest. Here, all Kji have an imaginary part. Equation ( 15. ip 
yields the fc-point correlation function by differentiation with respect to the source 
variables J in k [21 ESI El EH] , see below. The average over characteristic polynomials in 
general and their relation to determinantal structures are considered as well, see Refs. 

[SI S2l EH EU EH [Ml ES]. 

Let k 2 + h = k 2 , hi + 1\ — k\ and d = k 2 + N — k\ = Z 2 + N — l\ > 0. We consider 
probability densities P which factorize in the eigenvalue representation of the matrix 
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H. Due to the rotation invariance, we diagonalize H = UEU^ with a unitary matrix 
C/GU (AT). The measure is 

^ = M II 7-^ A ^)WM^), (5-2) 
where d[i(U) is the normalized Haar measure. Thus, we find from Eq. (14.171) 

^£L(«; A) = (-if^Nl n ^^zgjjK) (5.3) 

h/h j=l 

with 

sfo) = PiEjSfa) • (5.4) 
We decompose Zj into real and imaginary part, Zj = Ej + ijjj, and define the two sets 
k = diag («n, . . . , n kl i, K12, • • • , Kk 2 2) and A = diag • • • , n~ kil , Kk 2 +i,2, • • • , «fc 2 , 2 )- 

Our result for this choice, indeed, coincides with the one found by Borodin and Strahov 
|35j . They as well splitted the number of characteristic polynomials in two sets and 
derived the determinantal structure by discrete approximation. They used similar 
algebraic manipulations but they did not consider the connection to the supersymmetry. 
Hence our proof is truly a short-cut. As in Ref. [22], the splitting of fci and of k 2 in four 
positive integers is not unique. Thus, we find different determinantal expressions. 

We remark that we have only used the structure of the square roots of the 
Berezinians and no other property of superspaces. However, we may identify the terms 
1/(k o1 - n b2 ) in Eqs. fT4T27]) and (1428]) with the Efetov-Wegner terms [591 E3 E] which 
only appear in superspace. When calculating Eq. (15. ip with the supersymmetry method, 
such terms occur by a change of coordinates in superspace from Cartesian coordinates 
to eigenvalue-angle coordinates [60] [61] . 

The second term in Eqs. P~27P and (TOSjl . also contained in Eqs. @Og} and ( 14~29"j) . 
are intimately connected to the well known sum over products of orthogonal polynomials. 
This is borne out in the presence of M^ 1 which generates the bi-orthogonal polynomials 
[I]. Also, we might choose arbitrary polynomials in the square root of the Berezinian 
f!3.7p instead of the powers ■ If we take the orthogonal polynomials of the probability 
density P, then M^r becomes diagonal and Eq. (I4.26P is indeed the well known result. 
The fc-point correlation function can be derived by the case k\ = hi = k and l\ — l 2 — 
with k = diag {x\ + Liie — J±, . . . , + L k ie — J k , Xi + Liie + Ji, . . . ,x k + L k ie + J k ) 
where Jj are the source variables and Lj G {±1} with e determine on which side of the 
real axis the Kj are. The Cauchy integrals (I4.2ip and (I4.22p become integrals over Dirac 
distributions by summation over all terms with Lj = ±1 in the limit e \, 0. Thus, we 
find the orthogonal polynomials, too. Due to the differentiation with respect to the Jj 
at zero the Efetov-Wegner terms vanish and the well known result [1] remains. 

5.2. List of other matrix ensembles 

As we have seen for the ensemble over the Hermitian matrices, we find a determinantal 
structure (14 . 3 5 [) for particular matrix ensembles with help of the general integral ( I4.17P . 
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Here, we collect a variety of different matrix ensembles. Those ensembles share not 
more than two features: (i) the probability density function factorizes in functions 
of the individual eigenvalues and (ii) the non-factorizing part in the integrand is the 
squared Vandermonde determinant. We emphasize that this list is not complete. One 
can certainly find other applications. 

We introduce the decomposition into real and imaginary part, Zj = Xj+ti/j, and the 
polar coordinates Zj = rjC lipj . The probability density g(z) in Eq. (14.1 7p for particular 
ensembles with unitary rotation symmetry is up to constants listed in table [TJ 

Since the unitary ensembles describe physical systems with broken time reversal 
symmetry, one is also interested in ensembles which have orthogonal and unitary- 
symplectic rotation symmetry. For most of such ensembles the average over ratios 
of characteristic polynomials can not be transformed to one of the types of integrals 
discussed in Sec. HI However for the special orthogonal group and unitary-symplectic 
group and the Lie algebras thereof, they are integrals of the squared- Vandermonde type. 
They are listed in table El 

We remark that the integrals which have to be performed are different for every 
single matrix ensembles and can be quite difficultly to calculate. Nonetheless all 
averages over ratios of characteristic polynomials have the determinantal structure 
fl4.35p . The entries of the matrix in the determinant are more or less averages of two 
characteristic polynomials only. Thus, we achieve a drastic reduction from averages over 
a large number of characteristic polynomial ratios to averages over two characteristic 
polynomials for a broad class of random matrix ensembles. 

6. Applications for integrals of square root— Berezinian type 

In Sec. 16. 1[ we consider the Hermitian matrices again. We will show that the integral 
05. ip can also be understood as an integral of square root-Berezinian type. In Sec. 16.21 
we shift this ensemble by an external field Ho and transform the average in the usual way 
to an integral over a superspace. This integral is an integral of square root-Berezinian 
type. The determinantal and Pfaffian structure of the fc-point correlation function for 
intermediate ensembles is discused in Sec. 16 .31 

6.1. The Hermitian matrices revisited 

The integral (15.11) in eigenvalue-angle coordinates is invariant under permutation of 
the eigenvalues of the matrix H. As in Sec. [2J we present one of the Vandermonde 
determinants as a product over powers of the eigenvalues, 



N k 2 



N ( „ N 

3— 1 irajv i- 1 



n n (E a - km) 



,j-ln=l 6=1 





N fci 



n n (E a - K a ) 



a=l 6=1 



matrix ensemble 


probability density P 
for the matrices 


matrices in the 
characteristic 

1 * l 

polynomials 


probability density 


Hermit ian ensemble 

[53 EH E21 E21 El EE] 


P(trH m ,me N) 
Pf = #t 


H 


P(x)5(y) 


circular unitary ensemble 
(unitary group) 

jaTusnniaiiniEHiEniEiiEi] 


P(trP m ,m G N) 
C/tp = Ijv 


U and W 


P 


>^)5(r - 1) 


Hermit ian chiral (complex 
Laguerre) ensemble 
[65J EE E3 [7] 


P (tr(AAt) m ,m G N) 
A is a complex 
N x M matrix with N < M 


AA^ 


P(x) 




Gaussian elliptical ensemble 

m HUES]; for t = 1 
complex Ginibre ensemble 


exp 

x exp 
H is a 


(T + 1) trPttfj x 
2 J 

f (r - 1} RetrP 2 l 
2 J 

complex matrix; r > 


H and P f 


exp [— r 2 (sin 2 </j + r cos 2 </?)] 


Gaussian complex chiral 
ensemble [12] 


exp [-tr A ] A - tr P f P] 

C = iA + fiB 

D = iA^ + /xP 1 " 
A and P are complex N x M 
matrices with N < M 


CP and PtC 1 " 


Km-n ( 
x exp 


( 2^ rC ° S ^) 



cr 



o 
cr 

£ 2 to 



en 



matrix ensemble 



probability density P 
for the matrices 



matrices in the 
characteristic 
polynomials 



probability density 



O CD 

§ s, 

95 



real anti-symmetric matri- 
ces (Lie algebra of 
the orthogonal group) |1] 



P (tr P m , m G N) 
P = -P T = H* 
N = 2L + x dimensional 



H 



P(x)x x - 1 / 2 Q(x)5(y) 



03 P 

£ CD 

^ -a 

F! 3 



special orthogonal 
group [38j 



P (tr O m , m G N) 
T = 1 2L+X 



O 



P(x) 



--8{y) |1 — x\ x x 



x9(x- 1)9(1 -x) 



tr a- 

CO 
CB 



CD H & 



anti-selfdual matrices 
(Lie algebra of the 
unitary-symplectic group) 



P (tr P m , m G N) 



P 



P 



1 L 
-1l 



P 2 



1 L 
-1l 



P(a;)x 1 / 2 e(a;)5(|/) 



CD O 



w 

CD 



03 

£ t» 
CD 

* I 

™ CD 
- 



i2 

o' &■ 



unitary-symplectic group 
[381 



P(tr£ m ,ra G N) 



5 



1 L 
-1l 



5 



o m 

-1l 



P(x) Vl — x 2 5{y) x 
xQ(x - 1)0(1 - x) 
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Using the same decomposition of Zj = Ej + it/j as in Sec. El we identify this integral 
with the integral (14.1 1) and find 

3=1 V 7 

with 

f j (z j ) = Ej- 1 P(E j )S(y j ). (6.3) 

The integral (15. ip is permutation invariant with respect to the bosonic and fermionic 
entries of k. However, we do not see this symmetry in the expression found in Sec. 15.11 
because we split k into two parts. In the present section, we find a result which shows 
this symmetry from the beginning, see Eqs. (14. lip and (I4.12p for details. 

6.2. The Hermitian matrix ensemble in an external field 

Another calculation of integrals of the squared- Vandermonde type is not the only reason 
to consider integrals of the square root-Berezinian type. One of its powerful applications 
is the calculation of the A;-point correlation functions of ensembles in the presence of an 
external field. We generalize the result for arbitrary unitarily invariant ensembles over 
Hermitian matrices in Ref. [6] to ensembles in an external field. Assuming k < N, we 
consider the integral 



Herm(V) T ~ 



*■ det(H + aH — Kjit 



N 



where P is an arbitrary rotation invariant ensemble and Hq is an external field with a 
coupling constant a. For simplicity we set all imaginary parts of k equal to — e which 
means k = diag (x\ — ie — J\, . . . , — is — J/., x\ — is + J%, . . . , — ie + Jk) = x~ + J. 

We use the generalized Hubbard-Stratonovich transformation [61 [291 ED] to 
transform this integral to an integral over supermatrices. With help of the generalized 
Hubbard-Stratonovich transformation we arrive at 

Z = 2 2k{k - 1) J J $(p)Sdet~V + ® t N + al k+k ®H ) x 

x exp [— zStr p(a + + «;)] d[cr]d[p] (6.5) 

with cr + = a + islk+k- The superfunction $ is a rotation invariant supersymmetric 
extension of the characteristic function 



TP(K)= J P(H)exp(ttrHK)d[H]. (6.6) 

Hcrm (iV) 

Matrices in S (fc) are Hermitian and their entries are either of zero or first order 
in the Grassmann variables. The supermatrix p stems from the Wick-rotated set 
E^(/c) = U^T, (k)Il^ with the generalized Wick-rotation = diag (1^, e^' 2 lk) with 
ip e]0, 7r[. The supermatrix a lies in S_^,(/c) which is defined accordingly. The 
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Wick-rotation is for norm-dependend ensembles given by e*^ = i and guarantees the 
convergence of the integrals. For cases such as = exp(— Strp 4 ) we need the 

generalized Wick-rotation. Here, we assume that the probability density P yields a 
superfunction $ with a Wick-rotation regularizing the integral. 

We diagonalize the matrices p and a, p = UrU^ 1 and a = VsV~ l , and integrate 
over U and V which are in the supergroup U(k/k). Since we are interested in the 
/c-point correlation function R k of the shifted probability density P(H — aH ), we omit 
the Efetov-Wegner terms occurring from this diagonalization because they yield lower 
order correlation functions. The supergroup integrals are supersymmetric versions of 
the Itzykson-Zuber integral [HI [521 [61] . The integral (16.51) reads 



1 



(2m) 2k (k\ 



$(r)Sdet -\s + (g)t N + at k+k ® H 
BerJ%(«) 



det [exp(-w a is w )] 



Ka,b<k 



x det [exp(zr a2 s 62 )] 1 < a b < fc det [exp(-ir al (x b - Jb))]^^ 
x det [exp (ie^r a2 (x b + J b ))] 1<aM<k yjBer k %(s)d[s)d[r) . 



(6.7) 



l<a,6<fc 

Using the permutation invariance within the bosonic and fermionic eigenvalues of r and 
s, we find 



(2m] 



2 k 



<&(r) exp [— zStr rn] 



BergU^) 



x 



N 



nn 



e ^s a2 + is + aE i 



(0) 



exp [— zStrrs + ] y Ber k J k (s)d[s]d[r 



(6i 



a=l b=l 



Sai +is + aEl 

The integration over s is exactly an integral of the square root-Berezinian type (14. ip 



with the parameters N\ = N 2 
Fourier transform and find 



k, kx = and k 2 = N. In Appendix E we perform the 



/_ 1 U(fc-l)/2 



$(r) exp [— tStr rn] 



x det 




Bergi^) 



Ka,b<k 



X 



(6.9) 



oo 1 

T- 

^ n\ 



z?(o) 
iaEi V i 



n=N 



l<a<k 
Kb<N 



l<a<N 
Kb<k 



-aE, 



(0) 



a-l 



Ka,b<N 



S(r 2 ). 



We notice that the integration domain for the bosonic eigenvalues r a \ is the positive 
real axis where the integral for the fermionic eigenvalues is evaluated at zero. 

Indeed, we obtain the known result (HI EH] for non-shifted arbitrary unitarily 
rotation invariant ensembles for a — > 0. To show this, we put the 1/a terms of the 
Vandermonde determinant A^^aE^) in the last iV rows such that the lower right 
block is independent of a. The first iV terms of the power series of exponential function 
in the upper right block are missing. Hence, an expansion in k columns yields that up 
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to one term all other terms are at least of order a at the zero point. We find the limit 

(_l)fc(fc-i)/2 r r $( r )exp[-?Strr/d 
hm Z = - : / / , — - x 



(27T) 




x det 



al 



r a i - e^r h2 



5 



<9r 



b2 



N-i 



5(r 2 )d[r] . (6.10) 



l<a,fe<fc 

By differentiating the source variables in Eq. (16. 9p and setting them to zero, we obtain 
the modified fc-point correlation function 

k 



Rk(x~ 



K 2dJ j 




$(r) exp [— zStrrx" 
A iV (aE(°)) 



x 



(6.11; 



x det 



al 



r ol - e^r 62 



5 



2tt 



dr b2 



N-l 



oo 

ri l 



Ka,b<k I n! V 



l<a<fc 
!<6<iV 



z <9r 



l<a<JV 
Kb<k 



-aE. 



(o) 



a-l 



Ka,6<iV 



5(r 2 ). 



As discussed in Refs. [HI [29] , this correlation function is related to the fc-point correlation 
function over the flat Fourier transformation in x. Hence, we obtain 

r $(r) exp [— tStrrx] 



Rk(x) 



(2tt 



■2k 



xdet 



"2 



r„i - e^r 62 



A w (a£(fi)) 

JV-1 



x 



(6.12) 



-tip 



d 



2ty 



dr b2 



Ko,6<fe I * — ' n! 
< n=N 



oo 

2 i (^^ 0) ^i 



l<a<fc 
Kb<N 



% Or, 



1,2 



l<a<N 
Kb<k 



—aE, 



(0) 



a-l 



Ka,b<N 



We emphasize that this result is exact for any rotation invariant probability density as 
long as this integral above is existent. It generalizes known results [JTJ US] for norm- 
dependent ensembles. 



6. 3. Determinantal and Pfaffian structures for intermediate ensemble 



In Eq. ( 16.1 2p . we easily see that for factorizing characteristic function ( 16. 6 j) and, thus, 
for factorizing superfunction 



$(V) 



n 



$O i2 ) 



(6.13) 
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the /c-point correlation function is a ratio of a (k + N) x (k + N) determinant and a 
N x N determinant 



Rk(x) 



(2vr) 2fc A N (aEW) 



■det 



{i? 1 (x a) x 6 )} i ^ 



b<k 



R 2 {aE { b °\x a ) 



l<a<N 
Kb<k 



-aE. 



(0) 



a-1 



l<a<fc 
Kb<7V 



Ka,b<iV 



(6.14) 



Here, the entries are 



Ri(x a ,x b ) 



exp [- 



-i(rix a - e^r 2 x b 



N 



$(r)rf 

7™i — e % ^T 2 



-lip 



d_ 

<9r 5 



N-l 



/oo 
$(n) exp [-irircj — 



ri ) c?rj 



^a3(^b) 



2vr 



$ ( e ^r 2 ) 



exp [ze J ^r 2 Xfe] 



<9r<; 



a-1 



5(r 2 )dr 2 . 



6{r 2 )d[r] ,(6.15) 
(6.16) 
(6.17) 



By splitting off the lower right block from the determinant as in Eq. (12. 8p . we see 
that Rk is a A; x k determinant in x which was also shown in Ref. [41J. However, the 
representation (16.141) is much better suited for further calculations than for the k x k 
determinant representation. 

We can transform the characteristic function $ in Eqs. (I6.15p . (I6.16P and (16.171) to 
probability densities in an ordinary space by the inverse procedure performed in Sec. 16.21 
We emphasize that these correlation functions can also be expressed in terms of mean 
values over ratios of characteristic polynomials. To illustrate this we explicitly work 
out Eqs. (I6.15p . (16 . 1 6[) and (16.171) for Laguerre ensembles in Appendix F For Gaussian 
ensembles we obtain the correct result [46J. 

Instead of taking Hq as a constant external field, one can take it from a random 
matrix ensemble, too. For the Gaussian case this was discussed in Refs. l4"0t |4"TI |4"2"] . 
Here, we investigate Hq as a real symmetric, a Hermitian and a Hermitian selfdual 
matrix with factorizing probability density P. We remark that H can also be drawn 
from a Wishart ensemble since it can be mapped to one of the symmetric ensembles. 

Assuming that the characteristic function of P factorizes as well, the /c-point 
correlation function is 

i\JV(JV-l)/2 k 

■det[K^(x a ,x b )] 1<ab 



R ( k 2 \x) 



(2vr 



:2k 



,.b<k 



(6.18) 



for a second ensemble over the Hermitian matrices. We define the kernel 

K {2 \x a ,x h ) = R 1 (x a ,x b ) - 

RnsiXb) 



(6.19) 



N r ~ ~ 

/ P(E)R 2 (aE,x a ) (-E)™- 1 dE (M^- 1 ) 

1,72=1 TTj) 



ran qTi—1 
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Mi 2 2 



P(E)(-E) 



m+n—2 



dE 



(6.20) 



for probability density P. 

For quaternionic Ho, N = 2Q, we apply a generalization of de Bruijn's integral 



theorem |68J which we derive in Appendix C.2 This yields the Pfaffian structure 



4 4) (*) 



(27T) 



2k 



Pf 



~K^2 (xfe, x a ) -Kg (x a ,Xft) 



(6.21) 



l<a,b<fc 

where the sign of the Pfaffian for an arbitrary antisymmetric 2N x 2N matrix {D ab } is 
defined as 

1 n 

Pf [ £, afc]l<a,6<2JV = 2 ^ ^ II A^M^l) " ( 6 ' 22 ) 

w£6 2 jv 3=1 

The set Sm is the permutation group over M elements and the sign function yields 

and "—1" for even and odd permutations, respectively. The kernels in the Pfaffian are 

given by 

2Q 

K[ 4 \x a ,x b )= £ 



m,n=l 



a' 



mn Q,n—1 ' 



K?\x a ,x b ) = R 1 (x b ,x a ) + £ t} ^ L {M^) mn k^\x h 



2Q 



m,n=l 
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Ki 4 \x a ,x b ) = k^\x a ,x b )+ fcifw(^ ,_1 ) Ml *f ) ( 



x b ) 



(6.23) 
(6.24) 
(6.25) 
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The functions appearing in these definitions are 
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P(E) det 



-l) n / P(E) det 



9R 2 , „ N <9i?2 / r-i -v 

~^{aE,x b ) — (aE,x a ) 

R2(aE,x b ) E 
dR 2 



dE. 



n-l 



dE 



aE, x b ) (n — l)E 



n-2 



dE. 



(6.26) 



(6.27) 



In Eqs. f)6.23p . ( I6.24p and (16.251) the inverse of the skew-symmetric moment matrix 



M. 



(4) 



ah 



(b-a) P{E){-E) 



a+6-3 



dE 



(6.28) 



arises which generates the skew orthogonal polynomials of quaternion type. 

If H Q stems from an ensemble of (2Q + x) x (2<5 + x) rea l symmetric matrices, 
X G {0, 1}, we obtain another Pfaffian 



R^\x) 



f_l\N(N~l)/2 z k 



(27T) 



2k 



Pf 



K^\x a ,x b ) K^\x a ,X b ) 

-K^\x b ,x a ) K^\x a ,x b ) 



(6.29) 



Ka,b<k 
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The entries are 



2Q 



KV(x a ,x b ) = T (MW- 1 ) , 

1 v / j Oi m ~^ 'ran run— 1 ' 



m,n=l 



a' 



2Q 



Rm3 (^a) 



(6.30) 



(6.31) 



m,n=l 
2Q 



^(x.^) = fcW(x a ,o; b ) + £ feg)(x a ) (MW- 1 )^ (6.32) 

m,n=l 

with the moment matrix 

d[E] ,l<m,n<2Q + x 



-oo<Ei<_E 2 <oo 



-/ P(E)(-E) m ' l dE 



J P(E)(-E) n ~ l dE 







(-E 2 y 



\<m<2Q 
n = 2Q + 2 

X = 1 • (6.33) 

1 < n < 2Q 
m = 2Q + 2 

X=l 

m = n = 2Q + 2 



Here, the functions in Eqs. f)6.3ip and (16. 32ft are 

P(E) det 



and 



-oo<Ei<_E 2 <oo 



R 2 (aE 2 ,x a ) R 2 {aE 2} x b ) 
R 2 (aE 1 ,x a ) R^aEuXb) 



d[E] 



(6.34) 



<2(Q+x) 



/P(£)det 

-oo<_Bi<£;2<oo 



R 2 (aE 2 ,x b ) (-E 2 ) n - 1 
R 2 {aE x ,x h ) {-E^- 1 

-/ P(E)R 2 (aE,x b )dE 



d[E] 



l<n<2Q+ X 



(6.35) 



As in the other cases, the matrix {M mn } generates the skew orthogonal polynomials of 
real type with respect to P. Pandey and Mehta |39j constructed these polynomials for 
the Gaussian measure. They also found a Pfaffian structure for the interpolation between 
GUE and GOE. In Ref. |JT], one can implicitly recognize the determinantal and Pfaffian 
structure in the interpolation from an arbitrary Gaussian symmetric ensemble to GUE. 
Our results (|6.18|) . (I6.2ip and (I6.29P extend the determinantal and Pfaffian structures 
in the intermediate ensembles from an arbitrary symmetric ensemble factorizing in 
the probability density to an arbitrary unitarily invariant ensemble factorizing in the 
characteristic function. 
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Moreover, we can omit the factorization of the unitarily rotation invariant ensemble 
and find an integral representation in the superspace for an arbitrary unitarily invariant 
ensemble for an interpolation to the other classes of rotation invariance. For this purpose 
we integrate over H in Eq. ( I6.12p and find for the integral kernel a determinant or a 
Pfaffian determinant, depending on whether Hq is Hermitian, Hermitian self-dual or 
real symmetric. 

7. Remarks and conclusions 

We presented a new method to calculate mean values for ratios of characteristic 
polynomial in a wide class of matrix ensembles with unitary symmetry and factorizing 
probability density. This method is also applicable to the classical Lie groups and 
their algebras. Our approach is based on determinantal structures of Berezinians 
with arbitrary dimensions resulting from diagonalization of symmetric supermarices. 
Although we did not use any supersymmetry for ordinary matrix ensembles, we managed 
to reconstruct those Berezinians in the product of the characteristic polynomials with 
powers of the Vandermonde determinant. Using these determinantal structures, we 
obtained determinants whose entries are given in terms of the inverse of the moment 
matrix for the particular ensemble. These matrices are connected to the orthogonal 
polynomials and show that the known results from the method of orthogonal polynomials 
[23} ED ESI EU H] are obtained. In particular, we re-derived the results of Borodin and 
Strahov [32] for the ensembles over the symmetric spaces in a more direct way. 

The determinantal structure is stable when an external field is coupled to the 
random matrix. With help of the supersymmetry method, we have shown this for 
arbitrary unitarily invariant Hermitian matrix ensembles in an external field. Our 
formula for the fc-point correlation function is a generalization of recent results over 
arbitrary Hermitian matrix ensembles [5] and over norm-dependent ensembles in an 
external field [16]. Moreover, we considered an external field drawn from another 
symmetric ensemble. We calculated the /c-point correlation function for an interpolation 
between an arbitrary Hermitian ensemble factorizing in the characteristic function and 
an arbitrary symmetric ensemble factorizing in the probability density. We found 
determinantal and Pfaffian structures, too. For Gaussian ensembles, this coincides with 
known results [391 SI] • We gave explicit results for the Laguerre ensembles coupled to 
an external field in a way which is different from couplings investigated in Ref . p55J [70] . 
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Appendix A. Derivation of the Berezinians 



In Appendix A.l we derive the determinantal structure for (5 = 2. The derivations for 



the other two cases p < 2q and p > 2q for G {1,4} are given in Appendix A. 2 and in 



Appendix A. 3, respectively. 



Appendix A.l. The case (3 = 2 

Let p < q. The trick is to extend the ratio of products by additional 
variables diag (/Cp+i,ij . . . , Then, we apply Eq. (13.21) . Since these additional 

variables are artificially introduced and, hence, arbitrary, we perform the limit 
diag . . . , /t 9i i) — > 0. We find 

[] (K al - K bl ) Yl ( K a2 - K 62 ) 
l<a<b<p l<a<b<q 



V <? 
o=l b=l 

<? q 

n n ( K ai - ^62) 

a=p+l 6=1 



l<a<b<q 



P 1 

n (/«oi - «6i) n n (^i - 

p+l<a<b<q a=l b=p+l 

(_l)^+D/24« f[ ( Ka2 _ Kfel ) 
a=l fe=p+l 

n (^1 - /«6i) n n (^1 - 

p+l<a<b<q a=l 6=p+l 



EI ( K al - K b2 ) 
l<a,b<q 



det 



K al — ^62 



l^a,^? 



Kp+l,l = -" =K 9,l = 



det 



-1) 



q(q+l)/2+pq_ 



1 



K al — Kb2 



l<a,b<q 



n («oi - Kbi) 

p+l<a<b<q 



Sdet p - q K. 



(A.l) 



tp+l,l = --- = Kg,l=0 



This yields the result (13 .4p . 



Appendix A. 2. The case (3 G {1,4} with p < 2q 



Let p < 2q. This calculation is similar to the one in Appendix A.l The only difference is 
that we have to take Kramers' degeneracy in the fermionic eigenvalues diag (/«i2, • • • , ^2) 
into account. We first use non-degenerate entries to reduce the problem to the one for 
(3 = 2. Finally, we restore Kramers' degeneracy. We find 

n (K a l - K b l) EI ( K a2 - Kfe2) 4 
l<a<6<p l<a<b<q 



V 1 

II II ( K al - K b2) 
a=l 6=1 
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EI ( K al - «6l) II ( K a2 - «w) 
-l)/2 l< a < b <P l<a<6<2<j 



p 2q q 

n n (^ai - ««) n - K j+gt2 ) 

a=l 6=1 j=l 



^_ 1)g(g+ i)/2 +P A — 1 — det 



This directly leads to Eq. (13.111) . 



^P-2g 2g-p 
^al hb b2 

Kal - K b2 J !< a <P 
l<6<2g 

1 0/ J l<a<2g-p 
l<b<2g 



(A.2) 



«j2=tj+qr,2 



Appendix A. 3. The case (3 G {1,4} u>rf/i p > 2q 



Let p > 2g. Some modifications of the line of arguing in Appendix A.2 are necessary, 
we find 

n ( K ai - ft&o n ( k «2 - k&2) 4 

l<a<6<p l<a<b<i} 



p q 

II II ( K al - «62j 
o=l b=l 

n («ol - «6l) II ( K a2 - K fe2 ) 



(_-Qg(9-l)/2 !<«<fe<P 



l<a<6<2<? 



p 2q q 

n n - ««o n - «,-+^) 

a=l 6=1 j=l 



:-i) 



[p(p-l)+g(g-l)]/2 



X 



which implies formula (I3.12p . 



p-2g 
ol 



K>al ~ K b2 J !< a <P 1<«<P 
l<6<2g 1<6<P-2(J 



(A.3) 



'%2='%'+g,2 



Appendix B. Calculations of integrals of square root— Berezinian type 

Appendix B.l. The case k\ = fc 2 = k 



We calculate 



£N 1 +N 2 



i=i 



N 2 

n 



Ua^Uf^W^P^/N^mz] = (-i)^ n ^-^ x (B.i) 
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JVi N 2 



C N 1 +N 2 



3=1 



K al — ^62 J l<a,ft<fc ^ K a l — ^62 J l<a<fc 

l<b<N 2 



l<a<7V 2 -7Vi 
l<fe<fc 



<a<N 2 -Ni 
l<b<N 2 



1 



Zal — «62 J l<a<^Vi [ Z a l ~ -2&2 J l<a<iVi 
l<b<k l<b<N 2 



1 



We use the definitions (|4.4p to (14. 6p . The integral (IB. 1|) . then, reads 

iVi 7V 2 



y n 9j{zji) n 

, +jv 3=1 3=1 



(2) 

N 1 +k/N 2 +k 



[z)d[z] 



(_l)(^2+fc)(A r 2+fc-l)/2 det 



K al — ^62 J l<a,6<fc 

{G^/^)(K fe2 )}" i 



{F^)(« al )} 



Ka<fc 



M 



<b<fc 



N\/N 2 



(B.2) 



The next step is to extract the matrix MjvwjVa from the determinant. This yields 

Ni N 2 



3=1 



3=1 



= (-l) (7V2+fc)(JV2+fc - 1)/2 detM Wl/W2 det [i^^W,^ 
with K {Nl / N ^ as in definition (HTTjl . 

Appendix B.2. The case k\ < k% = k 
Let ki < k 2 . Then, we have 



Ka,b<k 



(B.3) 



-1 



l(fe 2 -fcl)iVi 



lim 



K*i + l,lvi K fc 9 ,l _ >' 00 



fc 2 

n 

j=fci+i 



N2-N1 7 (Ni/N 2 ) 
K n ^k 2 /k 2 1 



With help of Eq. (PP]) . we obtain 

7(^1/ 
J fci/fc 2 



,„ C_ 1 ^ 1 (fc 1 -l)/2+(fc 2 -fc 1 )7V 1 



k) (B.4) 



(B.5) 



x 



lim 



fc 2 

n 

3=fel+l 



det 



K 



N 2 -iVl+fc 2 -fcl 

3'1 



7 (ATi/Ar 2 ), > 
■^1/1 ^al) K 62; 

K al — K b2 



l<a,6<fe 2 



det [<r x ] 



l<a<fc 2 — k\ 
k L +l<b<k 2 



The limit expression is a function of 1/k i which is different iable at 1/k i = 0. Hence, 
using l'Hospital's rule we find Eq. ( 14. lip 
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Appendix B.3. The case k = k\ > &2 



For k\ > &2, we have to proceed in a similar way. We extend the number of the fermionic 
eigenvalues k<i and find 

ki 



Z $£ 2) ( K ) = {-iy kl - k *X N *- Nl \ liin J] k-.V"^Z: .7 - (/.') 



«fc 2 + l,2i"->Kfc 1 ,2— 



„Ni-N 27 .(Ni/N 2 )i 
v j2 ^fcx/Ai 



j=k2+l 



(_l)(te+2fci)(fc2-l)/2+(fci-fc 2 )(Af2-A f i) 



cf 5 1 7iV Ber t/*. ,/6) 



,(2) 



X 



X 



lim 



Kfei +l,lvAfc,,l -5-00 



fel 

n 



det 



jVi-jVa+fti-fca 

S'2 



1/1 1,^61, K 2j 



l<a,b<fci 



i=fc 2 +i 

This directly gives the result (14.121) . 



det [k^ 1 ] 



l<a<fci— k2 
fc 2 +l<6<fci 



(B.6) 



Appendix C. Extension of integration theorems for determinantal kernels 

Appendix C.l. Extension of Andreief's integral theorem 
We consider the integral 

{Sab} 



1= J det 

c N 



i r ab} 

l<a<k 
l<b<N+k 



l<a<N 
l<b<N+k 



det 



l<a<Z 
1<6<7V+/ 



l<a<N 
l<b<N+l 



d[z}. 



(C.l) 



The functions R a and S a are such that the integrals are convergent. Apart from this 
property they are arbitrary. We expand the first determinant in the first k rows and 
the second determinant in the first I rows and obtain 

V sign(p)sign(a)det[r ap(b) ] det[s aCT(6) ] 

' ' Ka,o<fc 1* 



k\(N- k)\l\(N-l)\ 



x 



a£& N+ l 

x det[R p{b) (z a , z* a )} det[S a{b) (z a , z*)} d[z\. 

J \<a<N \<a<N 

C N k+l<b<N+k l+l<b<N+l 

We apply Andreief's integration theorem for determinants [51] and obtain 

m 



Ka,b<l 



(C.2) 



k\(N- k)W.(N-l) 



- s ^{p)sign{a)det[r ap{b) ] det[s OCT(6) ] 

! -" ' l<a,b<k 1 



<a,b<l 



X 



x det 



R P (a)(z, z*)S a{b) (z, z*)d 2 z 



k+l<a<N+k 
l+l<b<N+l 



(C.3) 
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This expression is an expansion of a determinant of a (N + k + I) x (N + k + I) matrix 
in the first k columns and the first I rows. We find the final result 

{s ab } 



1= (-l) kl N\ det 



{ r ba} 



l<a<N+k 
Kb<k 



l<a<l 
l<b<N+l 

j R a {z,z*)S b {z,z*)d 2 z 

£ J l<a<Af+fe 

l<b<N+l 



For k = I = 0, we, indeed, obtain the original integral theorem by Andreief. 



(C.4) 



Appendix C.2. Extension of de Bruijn's integral theorem 
Consider the integral 

rr f A , \ {4*} {Ba(z b ,Z* b )} { C a ( Z b , Z* h ) } 

J = / det l<a<2Af+/ l<a<2N+l l<a<2N+l 

J L l<b<l l<b<N l<b<N 

C N 



d[z] . (C.5) 



As in Appendix CTT we expand the determinant in the first I columns and obtain 
1 ' 

" {S CT(o) (^,^*)} 



(2JV)! 
x / det 



x 



{C a{a) {z h ,z* b )} 

l+l<a<2N+l w /+l<a<27V+Z 

Kb<N Kb<N 



We define the quantity 

D a b = [ [B a (z, z*)C b (z, z*) - B b (z, z*)C a (z, z*)} d[z] 



d[z] 



Then, we apply the original version of de Bruijn's integral theorem [68J and find 



J 



(2N)\ 



sign(a)Y[A aU)j Pf [D a(aMb) ] 



l+l<a,b<2N+l 



(C.6) 



(C.7) 



(CA 



(TG&2N+1 

Summarizing all terms, the integral J is up to a constant 
{A ba } 



J-Pf 



{-A ab } 



l<a<2N+l 
l<b<l 



l<a<l 
l<b<2N+l 

Ka.b<2N+l 



(C.9) 



We fix the constant by the particular choice 

1/ 



l-A 



l<a<2N+l 
Kb<l 







(CIO) 



which yields 

J = (_i)tf(^-i)/2-H(!-i)/2jV!Pf 



{-A ab } 



l<a<2N+l 
Kb<l 



{A ba } 

l<a<l 
l<b<2N+l 

\<a,b<2N+l 



(C.ll) 



For I = 0, this is indeed de Bruijn's integral theorem. 
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Appendix D. Calculating integrals of squared— Vandermonde type 



We derive the cases (k 1 — k 2 ) = (h — h) < N and (/c 2 — &a), (h — h) < N in Appendix D.l 



and | Appendix D.2| respectively. 



Appendix D.l. The case (hi — fc 2 ) = (h — h) < N 

With help of Eq. (13. 7p . we rewrite the integrand ( I4.17P as a product of two determinants 



r N i 

/ n^w Ber £ 

4 i=i 



(2 /fe 2 +7v( 5 ) \/ Ber h/l 2 +N(z) d i Z ] 



"I) 



(h-fei)(il+fci-l)/2 



iV 



JJ^(z,-)det 



f 1 1 




r i i 






l<a<fci 




l<a<fci 



R2- 1 } 



l<b<fc 2 



l<a<d 
l<6<fc 2 



K" 1 } 



l<b<N 

l<a<d 
Kb<N 



x det 



{^} 



l<ft<Z 2 l<b<Af 



l<a<d 
l<b<l 2 



l<a<d 
Kb<N 



d[z] 



(D.l) 



Using the definitions (I4.20h - fl4.25j) . we apply the theorem of Appendix C.l and find 

N 



(2) , 
ki/ki+m 



z)\ Ber 



c N 



x det 



(2) 

h/h+N* 
1 



Xbl — A a2 J l<a<l 2 
l<b<h 



(z)d[z\ = (_i)(' 2 +^)ai+fc 1 -i)/2 iV! x 

{A d (A a2 )} 



l<a</ 2 



K al _ K 62 J l<a<fci 
l<b<k 2 



Zl/0 { K al, X 



1/0 
1/0 



l<a<fcj I J l<a<fc a 



{K d ( Kfe2 )} 



l<b</c 2 



l<fe<«! 

l<6<h 



(D.2) 



The last step is the same as in Appendix B.l We separate the matrix from the 
determinant by inverting it. This yields Eq. (I4.19p . 

Appendix D.2. The case (/c 2 — [h — h) < N 
We consider the integral 

N 

v i =1 



Ber ^ /ka+N (z) 



Ber 



(2) 

h/h+N 



{z)d[z\ = (_l)«i('i-l)/2+fel(*i-l)/2 x 
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N 



C N 



Y[g(zj) det 

3=1 



1 



K-bl — K-a2 J l<a<k 2 
l<6<fcl 

1 



Kbi — Z a J l<a<iV 
l<fe<fci 



Kf 1 } 



l<a<d K 
l<6<fci 



det 



1 



A51 — A a 2 J l<a<Z 2 
1<6<Zi 

1 



Aw - Z* J l<a<7V 

Mr 1 },. 



L<a<d A 
l<b<h 



d[z] . (D.3) 



Using the result of | Appendix C.l| we get 

N 



4 i=i 



(2) /fc 2 +iv( 2; )v Ber ! 1 2 // 2 +jv 



= (_l)('l+fel)(/l+Al-l)/2+AT(fc2+«2+l) x 



x N\ det 













A&l — A a 2 J l<a<Z 2 
l<ft<Zi 

{Am'} 



l<a<d A 
l<b<h 



1 



l<ft<fc 2 l<b<d K 



1/0 



l<a<fci 
l<b<h 



(DA) 



which is the desired formula. 



Appendix E. Calculation of the flat Fourier transform in Eq. ( 16. 8 j) 

We consider the flat Fourier transform 
J = 



a=l 6=1 



s i + is + aE. 



,-,„ ex P [-* Strrs+ ] v Ber fc/fc( s ) rf W • 



(E.l) 



By extending this integral with a Vandermonde determinant of — aE^ and using 
Eq. (13.61) . we find the determinant 



/_ 1 \fc(fe-l)/2 



Ui(r a6 )i _ {J 2 (r al ,a^ 0) )| 

Ka,o<fe I. J Ko 



{^3,a(r 62 )} 



l<a<7V 
K6<fc 



(0) 



a-1 



<fc 
KZK7V 



Ka,b<N 



(E.2) 



where s+ = si + le, s% = $2 + le^e and f 6 = diag (r a \, e^r^)- Hence, we have 
to calculate three types of integrals. The integrals in the off-diagonal blocks are the 
simpler ones. We have 

N 



J 2 (r al ,aE. 



(oh 



exp[-zr al s+] (-aE® 



(sf + aE- 



(oh 



dsi 



-i 



-aE, 



(0) 



N 



(N-l)\ 




t^exp 



-tr al sf + t(sf + aEf')ti + tsft 2 d[t, s 



Determinantal structures 



32 



) N+1 2n (-aE, 



,(0) 



N 



(N-iy. 

(-t) N+1 2ir (-aE-, 




S(h +t 2 -r a i)t$ 1 exp 



laE^U 



dt\dt 2 



N 



.(0) 

e(r«i) / if" 1 exp 



(iV-1)! 

- 27rz6(r a i) exp [«o;E 6 (0 V o i 

- 2m&(r al ) ]T - (laE^r^ 



iaEf ] (t 2 - r al 
1 

+ 2mQ(r al ) £ - (*a:£ 6 (0) r al ) = 



n=0 



i=AT 



and 



■hA r b2) = J (e # s^) a X exp [w fe2 s+] ds 2 = 2?r f- 

R 

The integrand of the integral 



-lf/> £ \ a 



dr b2 



-l 



£(r 62 ) 



Ji(r ab ) 



exp [-^Strr ab g + ] / 



A? 



d[s] 



(E.3) 



(E.4) 



(E.5) 



' s% — e~ l ^s 2 

has to be interpreted as a distribution. It is up to an Efetov-Wegner term 

J^fab) = 2vre^ J exp [-zStr f ab a + ] Sdet - N a + d[a] + 2m . (E.6) 

E-tf(l) 

This is the supersymmetric Ingham-Siegel integral [6j. We employ the result of 
Refs. (SI EH] and obtain 



Jl(fab) = -27r 



e(r al ; 



r ai - e^r b2 



—e 



-up 



d_ 

dr. 



N-l 



62 



5(r b2 ) . 



(E.7) 



Thus, we get for the integral ( 16 .ip 



J = det 



' al 



2tt 



d_ 

dr b2 
-l 



AT-l' 



l<a,fe<fe 



oo 

J2 h { l0cE b° )r al 



l<a<k 
Kb<N 



dr 



62 



l<a<N 
Kb<k 



-aE, 



(0) 



o-l' 



Ka,b<N 



X 



(E.8) 



(-2 7 r) fc (-l) fc ( fc - 1 )/ 2 0(r 1 ) ( 5(r 2 ) 
X A^pM) ' 

where 0(ri) indicates that every bosonic eigenvalue r a \ has to be positive definite. The 
distribution 5(r 2 ) is the product of all Dirac distributions S(r b2 ). 
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We consider the Laguerre ensemble 



N 



Pu(H) = n 



c\H c v+1 



exp(-ctrfi) det u HQ(H) 



(F.l) 



where i/, c G M + are some constants and is the Heaviside distribution for matrices, 
i.e. it is unity for positive definite matrices and zero else. The characteristic function is 

7P V {H) = c iN+u)N det- N - u (ct N - iH) (F.2) 

and the supersymmetric extension is, hence, 

$ v {p) = c^+^-^Sdet - N -"{cl kl+ha - ip) . (F.3) 

We notice that $ factorizes, i.e. it fulfills Eq. (I6.13p . Thus we can apply the calculations 
in Sec. 16.31 

The function R\{x a) x b ), see Eq. (I6.15p . is up to the Efetov-Wegner term, which is 
the normalization in this case, the same as the generating function (16.41) for k = 1 and 

a = 0, 

1 



Ri(x a ,x b ) 



x lim 

s\0 



X 



X a X}) 



(F.4) 



P»{H) 



det(H - x b l 



det(H -x b t 



N) 



Herm (TV) 

the 

'I 

Ri{x a ,x b ) 



det(H — (x a — te)l N ) det(H — (x a + te)l 



N) 



d[H]. 



Let the orthogonal polynomials of order N with respect to the probability density 
P u , i.e. tt^\x) = x N + . . . are the associated Laguerre polynomials. Then, we find 

7r( — I 



)N -l 2 



X 



X 



(N+p-iy. 

7r^ ) (cX )7T^i 1 (cX 6 ) - TT^^CXa)^ (cx b ) 



(F.5) 



cx a ) u exp(-cx a )Q(x a ) 



Xa X b 

which is indeed the determinantal kernel for the case a = 0. 

For calculating the second function R 2 (aEf\x a ), see Eq. (I6.16p . we consider the 
integral 



t MM" 

X n (aEl >,x a ) = i / $„(ri) exp -irizj : dr x 

J n\ 



It has a structure similar to Eq. flF.4j) . 



n(0) 



T n (aEf\x a )~ (aE ( b 

Pn (H) 



x lim 

e\0 



Herm (n+1) 



det(H — (x a — te)l n+ i) det(H — (x a + te)l n 



(F.6) 

(F.7) 
d[H}. 
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Thus, we have 

^°U>= (JV+ 2 ;! !),„, * (F.8) 

X ( tt B< 0) c)"ir;,';r _1_ " > (^)(CT a ) A ' + "- 1 -"exp(- Ca: je( a : a ), 

which also follows by directly integrating Eq. ( 1F.6[) . We combine this result with the 
definition ( I6.16P and find 

R 2 (aEf\x a )= (F.9) 
= [^ +V - X \c[x a - aEl 0) ])(c[x a - aE^])^ 1 exp(-c[x a - aEf ] ])Q{[x a - aEf ] })- 

2mc 



N 



(N + u-iy. 



- («^c) n ^r" 1 " n) (^ a )(cx a )^- 1 -"exp(-cx a )0(x a ) 

n=0 

We notice that the first term of R2(aEj ) \x a ) vanishes if x a is smaller than ctE^ . 

Also for the function R a 3(x b ), see Eq. f)6.17p . we find an expression of a form similar 
to Eq. tfEl) and Eq. f lR7j) . 

R a3 (x b ) ~ J P N+v+1 - a (H) det(H - x b t a ^)d[H] . (F.10) 

Herm (a— 1) 

We easily see that this is 

R a3 (x b ) = 2Txc 1 - a Tx[ N t u+1 - a \cx b ) . (F.ll) 

This result can also be obtained by performing the integration of Eq. (I6.17p . 

We emphasis that our result of the Laguerre ensemble in the presence of an external 
source is different to those in Refs. [63 [70] since the coupling is different. 
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